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LOCAL  GRAVITY  FIELD  MODELING 

INTRODUCTION 

Rffgicnal  or  local  grevity  fi«lci  models  arj  importent 
for  many  a?plica*ion s.  Csflacticns  of  the  vertical  and 
gravity  anomalies  are  required  for  alignment  and  correction 
of  nigh  accuracy  inertial  navication  systems.  Caflections 
are  re guirad  to  reduce  observed  horizontal  and  vertical 
angles  to  the  elllcsolo,  far  instance  in  apolying  the 
Laolace  correction  to  azimuths.  They  can  be  used  tc  convert 
between  astronomic  ana  j-sodetic  cooroinates.  Geoid  undula¬ 
tions  are  required  to  convert  betiaaen  elevations  aerivea 
from  satellite  oias  er  ve  ti  on  s  ,  for  examole  SFSt  and  heights 
above  sea  level.  undulations  crovice  a  direct  mro  cf  the 
geoid  in  the  region.  Thus  the  imocrtart  ouantities  to  be 
derived  from  a  model  are  gravity  anomaiy*  aelf>ctior3  of  dhe 
vertical  and  gecid  unauiations.  In  addition*  one  aculo  li'se 
to  be  able  to  usif  as  observations  in  concutir;  tn®  mcdel, 
all  available  geodetic  data  inclucin;;  as '.r  c~  ge  ode  t  ic 
deflections  of  the  vertical*  satellite  cerivea  -eoio  undula- 
tions*  observed  suavity  values  end  computed  snomsli??*  geo¬ 
detic  and  astronomic  coardinates*  and  soherical  h^  rTion  i  c 
geopotential  coefficients. 

Of  all  availaole  technioues*  only  least  siuarss  collo¬ 
cation  alloes  one  to  use  any  combination  of  r on-h om  ic aneous 
geodetic  data  and  to  ccxoute  any  of  th»  reiuirad  cuantitiss. 
In  addition*  least  scurres  collccation  gives  an  octimrl 
(minimum  variance)  methcd  of  i ne or pola ti ng  among  the  survey 
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data.  As  aith  any  pradicticn  nathad,  th j  rasolts  of  l»cst 
squarai  collocation  dapono  on  th-*  variation  of  th  a  srevity 
fiald  in  tn*  sraa  of  pridicticn  and  on  th«  dansity  and  accu¬ 
racy  of  tha  available  ob  s«rv  at  icns  .  in  rajiona  of  sHarply 
varying  topography t  a  »ajor  por*ion  of  tha  variation  of  th« 
gravity  fiald  is  dua  *o  this  topography.  Sinca  praoiction 
rosulta  ara  naarly  linaarly  rolatad  to  tha  "smoothntss"  of 
tha  ♦laid*  significant  inpreva-nent  5  can  ba  obtainsd  from 
laast  squaras  collocation  in  tnasa  casas  by  first  ramoving  a 
aortion  of  tnis  affact.  In  thj  cast*  this  has  oaan  accem- 
olishad  using  rathar  coarse  alavation  oodals  comoosad  of 
naan  alavations  or  elsvations  scalacf  from  topograohic  maos 
(“orsbarg  ano  Tscharning,  CPorsberg  and  •'adsan» 
19SC)#  (junkalt  193?).  fhs  data  spacing  cf  thssa  modals  has 
Jaan  on  tha  order  ot  Zj  :o  ZZ  arc  seccnos. 

This  reocrt  dascribas  an  axoariaant  conductad  to  ccm- 
outa  a  Iccal  gravity  field  nodal  using  laast  scuares  ccllo- 
cation  in  conjunction  with  savsral  digital  terrain  Todals  of 
varying  resolution  and  extant.  The  most  accurate  elevation 
data  teas  cbtaxnad  photogranmetrica  1  ly.  In  adoitionf  a  glo- 
oal  trenc  is  first  ramoveo  oy  computing  the  spnerical  har¬ 
monic  jxoension  of  tna  anomalous  potential  from  a  133  degree 
and  order  coefficient  sat*  .modified  by  the  kjSlZ  ellipsoidal 
ootantial. 
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BACKGROUND 


Th«  U.3«  Apty  Enginaar  fopograchic  La‘3orftori 
CU5AETL)  >^S3  ba«P  invalvad  jiith  davalacin?  larduara  and 
3oftJara  to  support  tha  jafa'^sa  «ap3ing  igarcy  CC'^A)  ana 
otnar  agancias  in  dansitying  ?oa  axtanaing  local  gravity 
nataorks  for  sany  yaars*  Vacantly*  this  aork  has  oeer 
diractad  primarily  to  aavalooing  inertial  survayinj  tachnel- 
ogy  for  datarmining  gravity  aromalies  and  cie'lections  of  t.na 
vartical  as  «a'll  as  positions  CTodd*  19?2>. 

This  aork  unit  aa s  initietad  aith  the  oojective  of  a*a- 
nining  othar  technologies  and  computational  methods  for  com¬ 
puting  magnitude  and  direction  of  th^  gravity  victor  and 
gaoid  unculeitions  at  a  large  rum  per  o*  points  in  6  local 
region  in  an  efficient  nay  and  aitm  sufticient  accuracy  for 
soma  or  any  of  the  aoplicatiohs  discusead  above.  In  cartic- 
jl3r»  computational  techniques  for  comeining  large  amounts 
af  heterogeneous  data  xera  ravieaad. 

The  interoolation  cf  gaodatic  cuantites  is  not  sc  much 
a  problem  in  flat  areas.  Tne  gravity  fields  in  these  areas 
are  smooth  enough  to  be  intersolated  adeauately  for  most 
ourposas  by  standard  1 oa s t-s quares  methods  Cdaiskanen  1  Mcr- 
itZf  1967).  The  da  term  in  a t ion  of  local  gravity  tield 
modalst  or  the  da la  nr  in  at lo n  of  the  gaoid  in  mountainous 
areas  is  still  very  much  e  prcoiem  hooever.  The  scarcity  of 
gravity  data  anc  the  great  irnaculerity  of  the  field  in 
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thas9  arses  make  tradition?!  iietnods  ina^fective.  ^or  these 
regionSf  least  squares  collocation  follosin?  trend  rjmoval 
oy  spherical  harmonic  expansions  and/or  digital  terrain 
node  Is  13  a  poitirful  tool  imich  has  haen  exoolited  in  the 
Nordic  nations  C=cr3oer5  *1  kasden*  19?C)t  Crorsherq  L 
Tscherning,  l)31a)t  C^orsderg  i  Tscharninj*  1951b),  es  a ell 
as  in  Europe  ( lurtner  L  rlmiger,  19S3),  (Sonkel,  19S3),  and 
in  Canada  CLachapelle,  1975),  In  addition,  )MA  has  con¬ 
ducted  in-house  efforts.  i  recort  on  these  efforts  CCarl- 
son,  1983)  was  obtained.  This  report  served  as  the  starting 
;)Oint  for  the  present  research. 

&  review  of  these  erforts  incicated  several  areas  in 
which  further  research  could  be  fccuseci.  These  include? 

1.  Effect  of  resolution  and  extent  of  dicital  t-»rr?in 
data  on  the  accuracy  of  the  solution. 

2.  I  nc  ro  ve  n.  n  ♦  5  to  be  gainec  by  using  a  very  dense  CT^ 
in  tre  iiraediate  vicinity  of  the  computation  point. 

3.  Effect  of  using  only  the  tf'*rain  and  isostatic  dste 
for  trend  removal,  neglecting  the  spherical  harmonic 
c  ompcnent . 

4.  Efficienf  methods  for  processing  the  terrain  data 
including  Fast  -ourier  Transforms  (=^T). 
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5.  f'athods  for  th*  inclusion  of  ^aoohysical  inforir?- 
tlcn  into  th«  procoss. 

THis  r*oort  sctnrassss  partially  items  Cl)  throuc:**  C  3) 
above.  Further  researcn  eill  oe  reouirec  to  more  fully 
answer  those  issues  ana  to  Tdoress  i*eas  C4)  and  (5)  end 
others  . 

The  next  section  briefly  describes  the  method  and  tne 
principles  oehina  it.  Tris  is  folloeed  by  a  description  of 
the  test  erea  and  data  used  ana  by  an  analysis  of  each  step 
of  the  process.  The  results  are  tnen  presanted*  follceed  by 
some  conclussions  and  reccmmenda ti ons. 
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DESCRIPTION  OF  THE  TECHNIQUES 


A  booy  at  rjst  an  th«  aartn's  supfaca  is  actad  upon  by 
tha  gravi tational  attraction  of  tha  aartb  and  tba  centrifu¬ 
gal  forca  causad  by  tba  aartn's  rotation.  Tna  effect  of  tha 
co.nb  in  at  ion  of  thasa  t«o  forces  is  callac  gravity.  Tha 
ootantial  of  gravity  In  is  than  the  sum  of  tha  gravitational 
ootantial  V  and  tha  centrifugal  ootantial  ^  . 

Thus: 


N  =  V 


1 


ui hare  t 


V  =  k 


a  no ; 


1  :  £  Z 

<t=  2avx  ♦  y) 


where; 


k  =  tha  gravitational  constant, 

V  *  the  volume  of  the  earth, 
f*  =  the  density  of  the  aartn, 

1  =  tha  distance  bataean  the  attracting  mass 
cno  the  attracted  mass, 
a  =  tha  earth's  rotational  velocity, 
x,y  *  tha  cartasi an  cooro*. nates  of  tha  point 
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und*r  consideration*  in  an  aartn^cant »r ad 
coordinate  frame. 


The  s'**vity  potential  can  also  'le  expressed  as  the  lum 
of  a  "nermal”  component  J*  enich  is  due  to  a  homogeneyaus 
ellipsoid  of  revolution  ehich  is  an  eouioo tent lal  surface" of 
a  "normal"  gravity  field  ana  a  "disturbing"  ootential  T, 
ehich  is  due  to  mass  anomalies  eitnin  the  actual  earth. 

Thus: 

e 

WCxtytZ)  »  TCK»y»c)  ♦  UCxtyti) 

In  Lj  IS  also  includea  the  cortrifugal  potential  and  the 
iffects  of  all  masses  external  to  the  eerth,  exterior  to 
the  earth*  v  satisfies  Laplace's  eauatien; 

e 

AV  -  0 

ehe  e  a  indicates  the  Laplacian  operator: 


Thus  it  is  a  harmonic  function  external  to  the  earth. 
Inside  the  earth,  it  is  not  harmonic  because  there  V 
satisfies  Poisson's  ecuation: 
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Th«  normal  potantial  U  can  oa  compiitad  diractly  for  a 
jivan  rafaranca  alliosoid.  All  that  is  raquirad  is  tha 
spacif ication  of  tha  four  quantitias: 

a»  samimajor  axis:  ft  flattaninj 

y.  t  aquatorial  sras/ity:  i'nd  a*  angular  velocity. 

Ha 

For  this  study*  tna  hGS72  rafaranca  allipsoid  aas  usad* 
tha  paranatars  takan  from  (CMA*  1974). 

Thara  ramains  than  tha  computation  of  T*  tha  "disturb¬ 
ing"  or  "anomalous"  potartial.  From  T  ua  can  dariva  all  of 
tha  important  guantitias  from  tha  aauations;  C-tetskanan  t 
•^oritif  1967) 


2 

gravity  anom?ly;  - -  "  ”  ^ 


aast-uia  st 

caflaction  componant;  a 


r  >  r  os  ^ 


d  T 

ix 


north-south 

Caflaction  comporant;  ( 


I 

ry 


dr 

d  <F 


gaoid  undulation;  N 


>a 


la 
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Th*  t«chniau«  untisr  i nv«s ti ;• t Ian  is  to  Dr*ak  T  into 
throo  parts.  Th*  first  is  a  ••global  aartn  nodal"  cowoonant 
dua  to  an  axpansion  of  tha  anon^lous  potantial  into  ?  sphar- 
ical  Harnonic  sarlas.  Tha  sacono  conponant  is  tha 
tooograp^lc-isost^tic  conponant*  computad  fron  a  ragienal  or 
local  nodal  of  tha  most  aall  known  mass  anomallast  tba  topo¬ 
graphic  massas  and  tnair  isostatic  conpansation.  A  third 
conponant  is  comoutad  from  a  sat  of  obsarvad  gravity 
anomalxas*  oaflactions  of  tha  vartical  and  gaoid  undulations 
in  tha  local  araa  of  intarast.  This  final  conponant  is  tha 
ona  actually  comoutad  ny  laast  squaras  collocation.  Thus 
tha  nodal  for  tha  anomalous  potsntial  is: 

T«TaT4T  11 
•M  t:  C 


iKhar a  ; 

^CM  *  T*!*  «srth  Tocal  conponant* 

T^  *  tha  topographic-isostatic  comoonant,  and 
Tq  *  tha  rasidual  cenponant*  or  daviation  from 
tha  modal. 

In  fraquancy  dcnain  tarns*  aa  may  consiJar  tha  "signal" 
T  to  ba  arokan  down  into  tha  Ion;  aavalangth  co'noonant  dua 
to  tha  spharical  harnonics*  a  nadium  and  short  aavalangth 
conponant  oua  to  tha  tcoograpy  and  its  conpansation  and  a 
rasidual  componant  dua  to  tho  daviation  of  tha  actual  grav¬ 
ity  fiald  from  tha  modal. 
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Th#p#  is  bowsvsr  a  ^roblan  a»ith  this  concaot  in  that 
ths  spharical  narmonic  coafficiants  contain  saraa  short 
aavalangth  information  as  jell.  shotld  therefore  remove 
the  t  opographic— “1  so  st  at  ic  affect  on  these  coefficients 
nefore  combining  them  with  the  other  components.  CForsherg 
i  Tschernlng  tl931a)  however  have  shown  this  effect  to  be 
negligible  for  fixea  area  computations  for  coefficients  up 
to  degree  and  order  36.  Morover*  the  effect  of  the  higher 
degree  and  order  coefficients  Jill  be  nearly  constant 
throughout  a  local  area  and  thus  this  effect  should  be  ccra- 
pansated  by  collocation. 

The  next  section  ciscusses  the  are»  itudied  and  the 
data  employed.  The  ccmputation  of  each  subcompcnent  of  T  is 
then  explained  in  detail  and  tn*  results  are  presentee. 
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TEST  AREA  AND  DATA 


Th«  crinclpal  t^st  lias  in  th»  aastarn  oart  of 
Vtvada  and  is  locatad  bataaan  $opro<isdta  north  latitudes 
37.5  and  3d.3  dagraas  and  asst  lon'^itudas  114.5  and  115 
dagraaa.  Tha  araa  is  sncan  xn  ♦igura  l.  As  can  ha  seen  in 
this  figurat  tha  araa  is  naarly  bisacted  ay  tha  Cry  Laka 
Vallay  and  is  fairly  Rountamous  outsida  of  this  valley, 
elevations  ranga  from  about  1400  matars  in  tha  vallay  to 
2000  natars  in  tha  mountains. 

From  aithin  this  araa*  a  sat  of  1344  gravity  obsarva- 
tions*  23  observed  astro- gaooetic  daflections  of  the  varti- 
cai  and  6  geoid  undulotions  obtained  at  doopler  observation 
iitas  aere  obtained  fron  CMA,  ill  of  tna  data  jiss  '•afar- 
ancad  to  WGS72.  -ron  the  :r?vity  observation s,  frae-eir 
gravity  anomalies  aere  cotainad  from  the  formula: 


AS  =g  *  .3C36H  -  jaema 
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ujhara : 


g  »  obsarveo  gravity  at  th?  »arth's  surf=tce 
.3086  3  frea-air  normal  gravity  gradient  in  mgal/meter 
h  =  heignt  of  th j  station  above  the  geoid 
gamma  =  ncrmel  gravity  et  the  corr asoond  in g  point 
on  tne  ellicsoid. 


A  subset  of  1062  of  the  ,ravity  anomalias  ahich  lie 
Jithin  the  principal  *cit  area  were  selected  for  initial 


li 


Figure  1 .  Test  area. 


analysis.  Tha  statistics  of  t*>is  data  ara  giwan  in  tabla  1. 


value 

hi  gh 

lout 

mean 

$ 

n 

Ag 

6  3.0 

-52.1 

-13.5 

22.70 

1C32 

ksi 

-1.^4 

-5.3 

-4.  i 

1.83 

23 

eta 

2.  54 

-12.4 

-4.r 

5.20 

23 

-24.69 

-27.05 

-25. 50 

1.C5 

6 

Table  1 .  Statistics  of  observation  data. 


Cigital  alavation  matPicas  covaring  tba  tast  araa  ana 
surrounding  araas  aara  also  obtairad.  Tha  coarsast  matrix 
has  a  spacing  of  1  osgraa  and  diaansions  of  10  dagraas  by  11 
dagreast  approx imat si y  csnter^c  an  tbj  principal  tast  araa. 
Th*  naxt  finsr  grid  has  spacing  5  "inutes  and  has  di mans ions 
3  by  6  dagraas*  again  ap  pr  ox  iita  t  »1  y  centaroc  on  tha  tast 
araa.  Tha  final  data  set  ccnsistad  of  alavation  values  for 
tha  princioal  tast  area  only.  These  elevations  xara  derived 
from  aarial  ohotograohy  of  the  tast  araa.  This  date  aas  rot 
given  on  a  grid  hut  aith  irregular  spacing.  ill  of  tha 
alavation  values  are  rafarancad  to  haan  Sa=i  Level.  The  1- 
degraa  and  S-mnuta  data  sets  «*ra  derived  from  topographic 
saps  and  have  an  astimateo  horicontal  sccuracy  of  100  to  300 
satars  and  vertical  accuracy  of  20  to  150  natars.  The  vert¬ 
ical  accuracy  of  tha  photo gramsatric  data  is  astimatad  to  oa 
under  10  faat#  oased  an  the  flying  height,  camera,  and 
staraoplottar  used. 


Tha  photogr ammatric  oata  aas  supplied  in  aopreximate 
aast-aast  profiles  aith  a  spacing  batuaan  profiles  of  530 


f««t.  Tha  spacing  bataasn  points  on  ths  profilts  is  not 
constant  Put  is  lass  or  «qual  tc  370  feat.  To  comouta  tha 
toooiraphic  and  isostatic  componants  of  tha  gravity  fiald«  a 
orogram*  dascrlbad  latari  aas  usad  anich  raquiras  tha  alava- 
tion  valuas  to  ba  givan  on  a  squara  grid  of  latituda  and 
longituda  valuas.  To  obtain  such  a  grid  fro.n  tha  ohotogram- 
satric  data«  softwara  aas  davalcpat  to  intarpclata  alavation 
valuas  at  tha  nodas  of  tha  grid  from  surrounding  valuas.  If 
a  data  point  aas  locatad  aithin  30  fast  of  a  noda?  tha 
alavation  of  tna  noda  uias  takan  to  ba  that  of  tha  data 
point.  Otharaisa*  a  laast  squaras  fit  of  tha  surrounding 
3lvatians  to  tha  aquation  of  a  Pilinaar  polynomial  in  x  and 
y  aas  parformad. 

To  astimata  tha  accuracy  of  this  intarpolatid  data  set* 
tho  alavations  supplied  aith  tha  1082  gravity  observations 
aa»*a  used.  These  aara  compared  mith  meightad  means  of  grid 
values  from  jiithln  a  carta  in  radius  of  tha  data  coints. 
oasad  on  a  comparison  of  106S  laall-snacad  valuasf  tha  RMS 
difference  was  plu3~sinu$  *3  feat. 

Tha  final  data  sat  usao  consisted  of  a  180  by  130 
degree  and  order  sat  of  fully  normalized  soharical  harmonic 
coefficients  of  tha  gaopotential  developaa  by  CPapo#  1981). 
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COMPUTATION  OF  THE  MODEL 


COMPUTATION  OF  THE  GLOBAL  EARTH  MODEL  COMPONENT 

Th0  •quatian  ♦or  r^prosonlim  th»  jartn's  gravity 
ootantial  W  oxtarnal  to  t^a  at'tracting  itassas  in  tarms  of 
sphtrical  harmonics  isJ  ("^aiskanan  C  *«oritz»  1957) 


kM  ao3nn_ 

y<rt^»A)  =  —  Cl  ♦  C-)  S  cosC.sX)  ♦  S  sinCmX))?  sin(^)3  ♦  lt>  (13) 

r  n»a.  r  m*c  nm  ns  nm 

ahara ; 

<t)  a  rotational  po+antial* 
r  =  jaocartric  dlstancat 
'i'  a  gaocantric  latitodst 
X  3  longituoa* 

kM  3  tha  jravitat icnal  constant  timas  t^a  mass 
of  tha  afrth* 

C  I  S  3  tha  fully  normalizad  soharical  harmonic  coafficiantsi 
nm  nm 

!*  3  tha  ♦ully  normalicad  sssocistad  Lagencre  p  o  I  ynotn  ia  1  s  « 

nm 

a  3  tha  acuatorial  radius  of  the  earth  and* 
n,m  3  tha  degraa  i-nd  order  r  esoactiv  ely  • 


In  oracticat  of  course*  equation  (13)  canrot  includa  an 
infinite  number  of  terms.  at  tha  orasant  time*  tha  oast 
solutions  incluoa  terms  up  to  ciegree  and  order  13).  Exam- 
olas  are  tha  Sodoard  earth  '<odel»  G:hl9C»  and  the  solutions 
of  Papp  (1973)  and  Rapp  (1951).  This  last  modal  was  used  in 
this  study.  To  obtain  T*  tha  anomalous  potential*  ee  must 
subtract  the  potential  of  the  reference  ellipsoid.  Tha 
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ootantial  of  tho  rofar^nco  •llipsoia»  'J»  may  ba  ojritten  by 


sotting  all  C  and  S  in  aquation  12  equal  to  zoro  except  the 

C,  ^  and  C,  ^  terns. 

2.0  4,0 

Thus: 


kM  a  4  34 

u(p,(^,x)  =  —  :i  ♦  c  — )  :*  p*  c  — )C'  p  i  ♦ 

r  r  ItO  it  0  r  4,0  *.0 

where! 

C'  *'/q  ■  potential  c oe t f ic iont s  conputad 

'  '  on  th»  basis  of  the  refererce 

e II  ipse  id. 


Subtracting  equation  Cl-*)  fron  eeuetxon  (li)  gives  the 
desired  Tg^  .  This  cjlculation  was  carried  out  fcr  all  of 
the  observation  ooints  usir;  tne  soherical  hamonic  coeffi- 
cunts  nodified  by  subtracting  the  i*GS72  elliysoioil  poten¬ 
tial.  The  latitude  of  tne  points  «as  first  converted  from 
geodetic  to  geocentric  to  conforn  to  eguations  (13)  ano 
The  values  of  a,  f,  k'*,  .and  C  for  tre  WGS 

ellipsoid  were  obtained  fron  COMA,  1974)  Tne  ouentities  of 
interest  ,Ag»{  » ’J  ,  and  M  ,  were  obtained  from  equcticns  T 
through  10  . 

For  the  purposes  of  tnis  study,  «e  are  oreeicting 
gravimetric  quantities  at  ooints  ohere  the  values  are  known. 
Thus  a  direct  measure  of  tne  accuracy  of  t^e  c  eirou  t?  ti  ons  is 
available.  For  the  case  of  computed  from  spherical 


(14) 
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narmonicst  th«  results  ar»  suO'iapi  zad  in  table  2. 


std.  dev. 


value 

?hS  orror 

m  ean  error 

of  error 

n 

A3  Cmgal) 

24. G3 

-6.40 

23.16 

1C82 

Ksi  (seel 

1.6e 

-C.2o 

1.62 

23 

eta  (sec) 

7.53 

-3. 24 

6.76 

23 

N  (m) 

1.27 

o 

• 

o 

0.31 

6 

Table  2.  Results  of  sperical  harmonic  computation. 


Given  the  snail  nunber  of  unci  ulat  ions  •  conclusions  are 
nearly  imoossible  for  t'^ese,  boeevar  the  solution  does 
closely  nodel  the  values  available  .  The  sravity  anomaly  vec¬ 
tor  is  not  well  nodeled  by  the  solution  hoaever.  One  oossi- 
ole  explar.ation  is  that  the  unculations  are  relatively  unaf- 
fecteo  by  the  short  javelength  tocobrachy,  ehtch  is  not 
nodeled  oy  the  solution  jhila  th a  jravity  vector  is  highly 
dependent  on  the  local  tooo3r?phy  ?nd  its  isostatic  compen¬ 
sation  . 


COMPUTATION  OF  THE  TOPOGRAFHIC/ISOSTATIC  COMPONENT  OF  T 

The  best  known  anc  oost  ?asily  observable  ness 
anomalies  are  those  associatjd  nith  th»  visiolj  toocgraphy 
of  the  earth.  Associated  aith  the  tooograony  is  the  "isos- 
i*^ic  compensation"*  or  tencency  of  the  tooographic  nasses 
to  be  compensated  by  mass  deficiencirs  oithin  the  earth's 
crust.  The  evidence  of  suco  ccopanscticn  is  given  oy  the 
oenaviop  of  deflections  ?nc  onjirelies.  esoecially  ioucer 
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anofflaliss  in  mountainous  ar«as.  Tha  Sou^a'*  rac'uction 
ramovas  tha  main  irragul ar it ia s  in  tJia  gravity  fiald  cssoci- 
atad  with  tha  visibia  toaography.  Thus  tha  vcluas  of  3ougar 
anomalies  should  ba  vary  small.  In  iioantainous  areas  how- 
avart  they  can  attain  values  of  several  hundred  milligals. 
4  similar  effect  can  be  observed  with  deflections  of  the 
vertical  calculated  solely  from  topographic  masses  in  moun¬ 
tainous  areas.  These  values  will  be  much  larger  than  their 
true  valuest  again  suggesting  a  comoensstion  beneath  tha 
mountains. 

ThuSf  not  only  the  visible  topography  but  the  gravita¬ 
tional  effect  of  its  isostatic  compensation  must  be  computed 
according  to  some  theory  of  isostacy.  T ra di ti on al ly t  the 
theory  of  Ai ry-H ei sk an er  has  been  used  and  this  theory  was 
^Iso  used  in  this  study. 

According  to  this  theory*  the  mountains,  j/ith  constant 
density/o  are  theorized  to  "float"  on  a  denser  layer  of  den¬ 
sity  P  .  As  the  mountains  are  in  floating  eju 11 librium, 
tha  higher  the  mountains,  the  deeper  ere  the  "roots". 
Analogously,  there  are  "anti-roots"  under  the  oceens.  The 
thickness  of  the  root,  t,  is  given  by; 

tAp  =  HP  CIS) 
■here  H  is  the  neignt  of  the  tooography  abcve  see-level 

and 
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is  dtnsity  contrast  bstpaan  tne  crust  and 


ths  Kant  la. 


3 

Usins  accaptad  naan  continental  v^luas  of  p  «  2.67  j/cm 

P 

and 

.  j 

Apa  -0.*  this  dacomas; 

6.C,7 

t  a  p  /O  K  H  =  X  h  (16) 

In  tf'is  stu3y»  jia  assume  a  crustal  thickness  of  32  <m 
for  tha  araa*  which  is  a;ain  consistant  with  tha  oast  kncwn 
continantal  values. 

Tna  terrain  haijhts  used  are  located  at  the  nodes  of 
regular  square  griast  as  explained  earlier.  It  is  natural 
then  to  compute  the  attraction  due  to  the  tooographic  masses 
and  isostatic  compensation  from  regular  rectangular  prisms. 
Tha  aidth  and  length  of  each  prism  is  simply  the  gnid  spac¬ 
ing  of  tha  particular  dtm  usao.  ^ or  the  topographic  effect* 
tha  haig.nt  of  tha  orism  is  the  height  above  sea-level.  For 
the  isostatic  effect*  the  height  is  the  thickness  of  tha 
"root"*  and  the  density  contrast  bataaen  tha  crust  and  man¬ 
tle  is  used  to  compute  the  attraction. 
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To  »nsur«  th*t  is  a  harmonic  function*  must  b* 
tha  potantial  of  a  givan  fixad  voluma  of  mass.  Thust  tha 
sama  voluma  of  tarrain  data  should  ba  used  for  each  cemputa- 
tion.  This  computation  icaally  should  axtand  to  such  a  dis- 
tanca  from  tha  computation  point  that  tha  tarrain  offsets  of 
tha  distant  topography  ramain  constant  for  th#  local  araa 
of  intarast. 

As  an  altarnativa  to  this  procadurat  ona  could  usa  a 
"rasidual  tarrain  modal**  or  tha  rasicual  batmaan  a  maan 
alavation  surfaca  and  tha  actual  tocography.  This  mould 
laad  to  significant  computational  savings*  as  shomn  by 
(Forsbarg  d  Tscharning*  1*315).  In  this  casa*  no  isostatic 
compan sation  naads  to  bs  cslculatad  end*  in  addition*  tha 
co/mputatich  can  ba  carrud  out  to  a  fixad  dist an ca  from  aach 
eomputation  point. 

Figura  2  shoms  a  reorasantati va  sjuara  prism  of  topo¬ 
graphic  cr  isostatic  mass.  Lat  this  prism  hava  constant 
dansity  p  ano  ba  attracting  a  point  mass  located  at  the 
origin.  Than  tha  vartical  component  of  the  attraction  of 
tha  prism  on  tha  point  mass  is  given  by;  C<allog*  1)53) 


.-'I: 


% 


y, 


vnara • 


G  =  tha  univarsal  gravitational  constant* 
r  a  distance  to  tha  attracted  point* 


(17) 


* 
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x#yf2  »  coordiratas  of  tha  diffarantial  altnant, 
P  >  dansity  of  tha  orisnif  anot 
r  >  potantial  of  tha  pris»> 


Intagration  of  this  formula  givas*  (Forsbarg  i  Tscharn- 
ing,  1981) 


i  a_rc.4-0Ln 


(S) 


N  y. 


CIB) 


Por  pristrs  ahicn  an  furthar  avay  ^rom  tha  attractao 
oointf  aa  nay  apcroxmata  tha  oriss  by  a  horizontal  ^sass 
olana  oassing  through  tha  cantar  of  tha  ppism  snd  oarallal 
to  tha  toe  and  bottcn  of  tha  prism.  In  this  casa?  the 
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sttpaction  in  tK*  vartical  direction  is  givan  by 


V, 


«nera : 

p'  *  eCz  -  r  ) 

1  2 

Z  s  <Z  >2  )/2 
mi: 

222  1/2 
rsCx  ♦y  +2) 


C19) 

(20) 


To  obtain  tn«  components  of  t*i«  attraction  in  the  hor¬ 
izontal  ciractionsf  th»  coordinate  exes  of  figure  2  ere 
rotated  about  the  origin.  Tnust  ee  use  the  seme  formules 
3ut  eith  the  t r ans f orm at ic ns ; 

t  t 

Cx  «x  ty  »y  tZ  iz  )  »  Cz  »z  »x  tx  .y  ,y  )  (21) 

121212  121212 

for  the  north-south  component,  endt 

t  t 

C*  *x  »y  .y  fZ  »z  )  s  (y  ,2  ,2  ,x  .x  )  (22) 

121212  1212  12 
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for  th«  9»sf-«9st  co^ipon^nt. 


Tn«n  th«  attractions  ar«  givan  by; 


3x.V  * 


C23) 


ahara  g*  has  batn  avaljstad  accord  in g  to  aquation 
C21)  or  C22). 


for  tna  conoansad  fcrmulast 
( 1 9>  givas  : 


intagration 


aouation 


C24> 


(25) 


Tha  Toco5r£phi.c  cr  isoststic  affact  on  th«  gaoici  uncu- 
lation  is  obrainjd  by  inragration  of  aquetion  (16)  mith 
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rtsptct  to  z 


N  =  (r  .( [»y  l-’a  ^ 


a.\rc.VaM/^\  (26) 


y. 


—  «Lre,ta.M  ^  -  -S.  a.fcta.n  I ~ 

Z  \y«-/  ^  \*r 


yi  *> 


x.v.*.' 


«or  distant  pns.sst  C  Pac '<i  1 1  an ,  1S5®)  ‘las  davalODad  simplar 
SKprossions  for  tha  potential  and  its  doriwatiwos  dy  a 
spnorical  harnonic  sxcansioo  of  th«  prlsi*  fiald.  Tha 
rasulting  har'sonics  sra  sinpla  polynomials  in  xy  yt  and  z» 
Tha  formula  for  tha  potantial  m  a  ccordinsta  systam  aith 
origin  at  t^a  cantar  of  tra  prism,  anc  •'xas  carp  jncticular  to 
tha  facas  is: 


y  *  s-asaxava* 


ZmT*  -  Ay*  -A*')x'  Ay*.  A2*  ^y* 


+ 


•h 


a«fr^ 


+  aaa 


C27) 


aharay  A*"»  Ayy  ^■'c  Az  sra  tna  sidalangths  of  tha  orisn. 
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An  «v»n  siinDljr  fcpro xinatzon  can  b«  vad«  far  vary  dis¬ 
tant  points.  Thasa  can  !}«  racnasantad  as  oelnt  ncssast  in 
ahich  casa  tha  pot-antial  is  jivan  ay  t^a  first  tarn  of  aqua¬ 
tion  (27). 

It  is  only  througn  tha  usa  of  tnasa  aoproxinata  formu¬ 
las  that  tha  conpytation  of  gravimetric  ouantitias  ♦rom 
digital  tarrain  nodals  bacomas  practical.  Soma  pralininary 
computations  vara  dona  using  only  tha  exact  fornulas  aith 
thraa  digital  tarrain  "odals  h?ving  incraasingly  aidar  spac¬ 
ing  amay  from  tha  computation  ooint.  Thasa  aara  tha  ona 
dagraa  anc  fiva  minuta  matrica>  dascribad  aarliar  and  an 
additional  thraa  sacond  naan  al.vation  matrix.  A  total  of 
46#494  prisms  must  ba  awaluatad  for  ooth  tha  topographic  and 
isostatic  comoonants  ahan  using  tha  ona  dagraa  maan  alava- 
tlon  matriXf  tha  fiva  ’iinuta  ra?n  alavation  natrix  ano  a  tan 
minuta  by  tan  minuta  ragion  of  thraa  second  maan  alvetions. 
This  computation  toot<  fporoxitr  at  »ly  ona  hour  oar  point  on 
tha  VAX  11/73C  comoutar  systamy  inclucJing  cpu  tima  and  I/C 
tima.  Tha  sama  computation  took  eoproxinataly  rf  minutas  par 
point  lihan  a  combination  axact  and  approximsta  fcrmulss 
a art  usad.  Tha  axact  prism  fcrmulns  wara  usad  cut  to  a  dis- 
tanea  of  ^  o  from  tha  computation  point*  ahara  d  is  tha 
langth  of  tha  main  diagonal  of  tn#  prism  baing  usad.  Tna 
'AacHillan  formulas  aara  than  usad  out  to  a  aiitanca  of  d 
from  tha  computation  point.  At  cistancas  further  than  this* 
tha  point  mass  ^crmjlo  j*a3  usao.  Tnasa  cut-off  distances 
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rtav*  b»tn  datarninad  by  1944)  ta  .iva  e  goad 
trad«-off  bataaon  tinia  aoe  accuracy  and  aara  varifiad  during 
this  axparimant.  An  axe  action  ta  tnis  is  that  tha  a ass 
olana  fornulas  aara  alaays  usad  for  gaaid  undulations. 

Tna  topofiraphxe  and  isostatic  affacts  aara  first  :alcu- 
iatad  for  a  subsat  of  tha  odsarwations  using  a  variaty  o* 
alavation  oata  and  taenniouas  to  gat  a  faal  for  tna  traoa- 
offs  in  tinia  and  accuracy  mvoLvad.  In  aadition  to  tha  23 
ooints  ahara  daflections  aara  knovn*  52  points  anara  gravity 
has  baan  obsarvad  aara  us?d  ss  tast  points.  Tha  ana  dagraa 
naan«  fiva  minuta  naant  ?na  fiva  sacond  point  alavation 
Tatneas  aara  usad  alona  and  in  cosbination  xith  aach  other. 
Tha  topogr  apni  c/isos  t?  ti  c  componants  j»ara  than  subtracted 
from  tha  fraa-air  anomaiias  and  kno^n  daflactions.  Tna 
results  sill  assantislly  ba  tha  "isostatic  anowalias  " 
Chaiskanan  i  *'oritz»  IS57)  in  the  Cusa  of  anoaelias  end 
“topographic  isostatic  aa  f  lac  t  ion  s '*  C'^siskanar  1  Mainesz)  in 
the  case  cf  oaflectioni.  Tha  results  are  oresanted  ir  t?bla 
-.  Tha  resulting  anorcalias  ond  dafl*ctions  should  theoreti- 
:ally  b»  small  end  have  lo«  vsrifness.  ’’sola  2  can  than  be 
used  to  judge  the  comparative  at  feet iv anes s  of  using  various 
terrain  date  sets  by  coirpering  the  maans  ano  standard  devia¬ 
tions  in  this  taole. 

Several  interesting  observations  can  oa  n ?c j  frot  this 
tabla.  Tha  inclusion  o*  the  fiva  second  data  -arginally 
improvad  tna  results  m  this  cass.  hora  imcortant  than  *his 
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Terrain  3ata 
Jsed 


4  3 


Value 

As 

mean 

5.73 

-5.5  5 

high 

Z^,Z2 

-1.33 

loe 

-7.^3 

-12.41 

a 

7.30 

3.03 

mean 

-0.14 

-1.87 

high 

j.  73 

0.5  ’ 

loa 

-4.71 

-4.56 

(7 

2.05 

1.20 

mean 

-3.83 

-0.07 

high 

3.03 

4.3  3 

loe 

~  1  i  •  :  6 

-6.52 

ff 

5.1C 

3, 4  < 

0 


Cl  ) 

Cl) 

5. 33 

0.20 

-5.54 

-0.85 

24.13 

4.44 

-0.58 

4.09 

•5.66 

-6.39 

-12.25 

-7.62 

5.65 

3.04 

2.97 

2.99 

C2) 

C2) 

-1.30 

-1.79 

-1.57 

-1.63 

0.18 

0.66 

0.40 

C.49 

-5.00 

-4.46 

-2.20 

-4.76 

1.18 

1.20 

1.01 

1.20 

C2) 

C2) 

0.11 

-0.11 

-0.11 

-1.62 

4.97 

4.34 

4.33 

2.20 

-’.16 

-6.55 

-6.52 

-6.2  2 

3.7? 

3.48 

3.51 

2.63 

Table  3.  Observed-component  from  topographic/isostatic  data. 


Notes: 

(1)  .iased  on  n  =  4Ci  as  12  points  fell  outside  the  innermost 
terrain  data  grid, 

C2)  dased  on  n  =  22,  as  1  ooint  fall  outside  the  innermost 
terrain  data  grid. 

Terrain  data  grids  used: 

4  -  One  degree  mean  elevation  grid  only. 

3  -  Five  minute  mean  elevation  grid  only. 

C  -  rive  minute  mean  elevation  grid  and  no 
bicubic  spline  interpolation. 

3  -  One  degree  and  five  minute  elevation  grids. 

5  -  Five  minute  mean  and  five  second  mean  alvation 

grids. 

^  -  One  degree,  five  minute  and  five  second  grids. 
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innar 

grid 

aas  tha 

aoi li ty  to 

Construct  a 

vary  dense 

"innar 

-inr ar 

"  grid  by 

interpolation 

of  tna  given 

elevations 

in  a 

three 

by  three 

grid  soacing 

ragion  of  the 

data  point. 

This  «as  4econpli»n»d  by  a  bicuoic  sclln*  int >r so latibn  of 
thosa  alavations.  This  intarpolsf ion  jas  usad  for  ?11  casas 
axcaot  for  ona  casa  in  ahicn  only  tha  fiva  ainuta  maan 
alavations  jiart  jsad*  Tna  iaprova'aant  can  ba  sa^n  by  com- 
oaring  tha  rasults  obtainad  from  alavation  data  sata  3  and 
\ot#  that  tha  rasulto  ara  imorowad  for  gravity  anomalias 
out  not  for  daflactions*  carnacs  r^jflactinj  that  daflactions 
art  lass  affactad  by  tha  wary  local  tooograpny.  Anothar 
intarasttng  obaarvation  is  tHat«  again  for  anomaliaSi  tha 
jsa  of  tna  ora  eagraa  oaon  aiavations  halos  to  bring  tha 
naan  rasiCual  vary  naar  to  2sro»  This  carhaos  indicatas  tha 
dasirability  of  using  such  s  broad  and  coarsely  soacad  grid 
for  th»  outer  zones  tc  rsmovj  local  oiasas  incurraa  from  tha 
innar  zona  oat.’. 

For  tnasa  cas?«f  tha  most  Casir-ble  conbinaticns  of 
alavation  data  socaar  to  oe  tna  one  CJacr>a  ^rid  aith 

tna  five  irinuta  narn  jridt  using  an  intarpol  at  ion  to  dansify 

tha  fiva  minuta  grid  in  tha  vicinity  of  tha  comoutation 

ooint  and  this  satra  date  set  uitn  also  the  *iva  saconc 

siavation  dstat  or  cate  sats  0  and  occor&ing  to  tabia  3. 

To  orovide  a  wida  variety  of  casas  for  ♦ha  laest 
squares  collocation  program,  another  larger  sat  o^  observa¬ 
tions  uara  used  as  test  points.  Using  the  ona  negraa  and 
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at  500  points  Jithin  tast  ration.  Csin;  tnis  alavatlop 
data  togathar  Jith  th^  fiva  sacond  point  aia^atiop  grid* 
anomalias  aara  cojaputad  for  a  suosat  of  110  of  t^asa  poin-ti. 
?oth  sats  of  tast  points  aara  c'losan  to  ba  i/all  distributad 
spatially  and  to  hava  tha  sama  alavation  oistogra'n  as  t^a 
original  ICo2  oolnts*  Th>  avaraga  sp^cin;  batsaan  tha 
points  aas  anout  1  arc  minut*.  In  ahdition«  tna  undulation 
aas  conputad  at  tha  sIk  cocol^r  points  usinj  tha  sacond 
alavation  data  sat.  Tha  locations  of  tha  tost  points  ara 
shoan  in  figura  3. 

Tha  valuas  co.sputaa  troa  tha  tarrain  data  sera  than 
subtractad  from  tha  fra*  sir  gravity  anomslias  and  obsarvad 
undulations.  Tha  rasults  are  shown  in  table  4, 

TerrsmlataL'ssd  0  “ 


n 

'  )C 

110 

Te  .'r 

-11.  ‘2 

-14.32 

;n 

: 

0.43 

1  :tt( 

-3*. 33 

-3^.50 

1  = .  3  7 

17.97 

?.  5  7 

10.50 

Table  4.  Gravity  anomaly  residuals  (observed-computed  from  topographic/isostatic  data. 

Terrain  Cats; 

I  -  Cna  d»gr«j  jnci  fiv?  rinjts  se  an  alevations. 
r  -  Cn»  Ci«gres«  ♦iva  -ninutei  and  fiva  sacond  elvation 


For  geold  undulations,  the  one  degree  and  five  minute 


30 


2290 


2290 


2270- 


Lat 


2260- 


2250- 


I  I  I  I 

14700  1471C  14720  14730 

Lon  (both  lot.  and  Ion.  tro  in  ninutts> 

figuro  3.  Observation  data. 

Logond: 

.  a  Crowity  Anomaly  Oboorwation 
a  a  Co71octien  Obaorvatlon 
e  a  Cocolon  Station 
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data  iats  (T«r'*3in  :ata  -i  o'  taole  3)  aara  ji«d.  Tha 
rasuits  jara  poarp  aitr*  a  Tiaan  «rror  of  -^.03  natars  ^nd  a 
standard  daviation  of  tna  irror  st  around  0.9  'latars.  A 
oossibla  axplanat-ion  is  t!ia  "indiract  affact”  of  ti^a 
topograpnic/isostatic  raduction.  H»).si<anon  and  '*oritz 
C1947)  stata  that  this  atfact  can  la  on  tha  ordar  of  IC 
aatars.  It  aould  ba  s  rslativaly  simpls  aattar  to  chack  on 
this  by  computing  tha  potantisl  dua  to  tha  tooojraphy 
saparataly  from  that  computad  from  tha  isostatic  comoansa- 
tion.  Cividing  tha  aiffarjnca  batu»a»n  tha  t  j»o  by  normal 
gravity  aill  giva  tha  indiract  affact  }n  geoio  undulation. 
This  aas  not  dona  in  this  study. 

COMPUTATION  OF  THE  RESIDUAL  COMPONENT 

Tha  rasidual  conponantt  ■*"*«  component  net  modalad  Oy 
aithar  tha  glooal  aartn  acoel  or  tha  topograohy-isostecy  a: as 
conputad  by  least  siusros  collocation.  Tna  tooographic- 
isostatic  comoonants  of  gravity  anon alias#  dalfoctions#  and 
undulations  ivara  ac/dsd  to  tna  corrasoono  ir  j  glcoal  earth 
nodal  comoonants.  Thjsa  suns  uiara  than  subtracted  from  tna 
obsarvad  cata  to  give  "rasiduel"  cosarvations.  Th.se  uera 
used  to  compute  rasicual  values  at  other  known  ooints  which 
jiara  not  usad  as  oos erv ati on s.  The  predictad  residuals  vara 
then  comparad  with  tha  knovn  resiauals.  ^he  statistics  of  the 
residual  observations  are  given  in  table  5. 
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Tarrain  Data 


0 


V  alua 


n 

soo 

lie 

an 

-4.74 

•"7.3  5 

f'i  -:b 

14.35 

9.10 

1  0  « 

-25.0  7 

-27.63 

.2 

■>2.  9  0 

33.71 

n 

2  3 

23 

■na  an 

3.2  1 

3.47 

ni  ,h 

5.  36 

5.89 

lev 

0.3  7 

0.57 

a2 

1.5  3 

1.76 

2  3 

23 

re  ?r 

C  .  9  2 

-0.61 

hi 

5.23 

3.14 

1  0  V 

5  7 

-i  .23 

11.  7. 

7.0  3 

n 

o 

Tie  »n 

-8.0  3 

(not  cjmoutad) 

bi  ;h 

-7.13 

1 0  a 

-9. <.2 

0.6  0 

Table  5.  Residual  values  (observed-(TI  &  EM)). 
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Th«  soal  of  raducinq  tha  ob  st  r  vat  lo  s  oa'ora  acol/ing 
collocation  was  to  snooth  out  tha  gravity  fiald«  or  raduca 
tha  rang*  of  valuas  and  th#  varianca  of  tha  gravity 
anomaliasa  Froia  -  tabla  5t  aa  saa  that  tha  varianca  of  tha 
anomalias  has  baan  raducac  fron  470.96  to  oataaan  ^0  and  )Q 
daoanding  on  tha  alavation  data  usad.  Comoarison  jiith  tabla 
1  shoas  that  tha  varianca  of  all  ot  tha  ouantitias  has  oaan 
raducad.  Anothar  dasirabla  attributa  of  a  "snteotn"  gravity 
fiald  is  that  tha  naan  valuas  of  tha  ouantitias 
g*  (  »  >7  «  and  N  oa  naar  to  2aro.  Tha  aaan  of  tha  unra- 
ducad  anoaalias  aas  -14.60  itgals  ahila  that  of  tha  raducad 
valuas  is  -4.74  agal  for  tha  bast  casa.  Tha  naans  of  tha 
daflaction  componants  ara  also  brougnt  naarar  to  zaro  by  tha 
raduction  and  so  is  that  of  tha  undulations. 

Thasa  rasidual  obsarvations  aara  than  usad  to  form  the 
observation  eguations  fcr  lasst  souaras  collocation.  Tha 
oasic  form  of  an  observation  eguation  uitnout  oarameters  ist 
(horitz.  1972); 

K=s»n  (2i) 

a  nara : 

X  *  a  vac  tor  of  oosurvations*  in  this 
casot  residual  gravity  anomaliest 
deflections*  and  geoid  unculations* 

s  a  tne  unknown  ••signal"  oortion  of  tha 


oas4rv«tiont  sndf 


n  *  a  random  "noisa"  comeonant. 


Tha  euantitias  s  ano  n  ara  assuaad  to  b«  unc orralatad* 
zaro-maan  quintitiaa  aith  Gausaian  distributions.  To 
oradict  tbs  signali  Tt  at  otner  unknoan  paints*  tha  Iszst 
squaras  coileeatian  foriruls  for  tha  minimum  varianca  asti- 
mata  is* 


*-  -- 1 

s  »  C  C  X  C2'9) 


ah  ar  a ; 

«» 

C  s  cov(s*$)f  tn»  matrix  of  covariances 
s  betaaan  rha  quantity  being  estiiaatad 
£.nj  tha  ab sarvat ions »  and 

C  =  CcovCsts)  ♦>  ;ov(n»n))*  the  sum  cf  tha 

matrix  of  ccwariancas  amang  and  bataaan 
♦he  qh servrt ions  end  tha  cavarianca 
•"etrix  cf  the  noisa  vector. 


Sinca  the  quantities  ae  jish  ta  estimate  ere  linear 
functionals  of  tha  anomalous  patential  T,  it  suffices  to 
datarmina  an  expression  for  the  spatial  covariance  of  T  and 
then  to  apply  tha  laa  of  preoojation  of  covariances  Choritz* 
1S72)  to  this  expression  ta  obtain  the  covariances  betsaen 
all  of  tha  observations  and  predictions.  Let  KCPtC)  be  the 
covariance  batmaan  T  at  point  ^  and  T  at  point  C.  Noa  lat  a 
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3nd  t)  b*  the  t«3  quartitias  darivec  froi  1*  by  linaar  scara- 


tlans  sue^  as  dif f ar an ti at  ion  or  multiplication  by  a  con¬ 
stant.  If  ■«  Mish  to  symboliza  thasa  oparations  by  KT)  and 
dCT)»  than  tha  cov-arianca  oataaan  a  and  b  is  sivan  by: 

r  Q 

CCa  lb  )  «  covCa  ,0  >  «  i(3(K(P,l}))  C30) 

P  I  0  C 


► 


Tha  laast  squaras  collocation  proqram  of  CTscharninqt 
1974)  aas  usad.  This  program  allouis  tha  sc  ac  if  Ic  at  ion  of 
ona  of  thraa  isotropic  (azimuth  indaoandant)  covarianca 
functions  ahich  ar a  all  basad  on  on«  of  tha  anomaly  dagraa 
sodals  davalopad  by  (Tscharning  £  5aoo,  1974).  Tho  "dagraa 
variancas"  of  tha  jravity  anomalias  ara  givan  by  C^atskannan 
and  Moritz*  1967); 


c 


n 


n,v 


m  "O 


nn 


whara ; 


(31) 


a  ffc  =  ♦ully  nornflizad  spharical 

nm  nm  her-^onic  cc*fficiant«  of  tha 

5'"awity  anosalias. 


Wa  can  than  arita  an  axprossion  for  t^e  covarianca 
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tunctian  or  tha  jravity  aromaiias  as; 


CCPfJ)  »  cavCAg  ,^3  )  » 

«  C 

iah«r« : 

P  s  th«  l»j»ndra  polynomial  of  dagraa  n* 

^  3  tha  spharical  distanc*  bataaan  P  and  Of 

R  3  radius  of  tha  Bjarnamnar  SDhara« 
r  3  distanca  from  P  to  tna  origin*  and* 
r*  =  distanca  from  0  to  tba  origin. 


R 

y^  C  P  Ceos  i/.  )C - 


n»0 


n  n 


R  n*2 

> 

rr  '■ 


Tho  dagraa  variancas  of  tha  anomalous  potantlal  art 
ralatad  to  thosa  of  tha  gravity  anomalias  by; 


l<  3  - 2  c  (33) 

n  Cn-l)^  n 

wnara; 

R  3  tna  maan  aarth  radius. 


NO  can  arita  an  anprassion  for  tha  covarianca  of  tha 
anomalous  potantial  as: 


<C?*0)  =  t( ^  ) 


Z 

ns© 


K  P  Ceos 
n  n 


) 


C34) 


Tha  modal  usad  for  tha  anoaaly  dagraa  vsrisneas  aas; 


I 
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C32) 


a  (1-1) 

c  «  -  (35) 

1  (i-2)(l^e) 


«h*r«  th*  constants  A  and  3  ara  datarmin^o  fret 
tbo  varianco  of  tho  obsorvad  cravity  anosaIi)s. 


Tho  itothod  of  collocation  olso  oivas  an  •stivata  of  th* 
naan  squara  arror  of  praoictxon  asi 


2 

«  3 

s 

«hara: 


t  _-l 

c  -  C  C  C  (36) 

S  S  }  s 

C  3  cov(»  »s  )«  tha  oolnt  varianca 

S3  ?  ? 

of  tha  quantity  jsti.iatad. 


In  this  casa«  hoaavar,  ««  kncu  anactly  tha  arror  of  tba 
astimata  sinca  aa  ara  cracicting  at  knonn  ooints. 

Savaral  casaSf  irvolwir”;  diffarjot  suosats  of  obsarva- 
tions  ana  pradictionSf  aara  run,  Tha  observation  ooints 
«ara  saloctad  to  have  on  avaraja  soacinj  cf  scop  ox  iirat  aly  5 
arc  ninutas.  Tha  oate  used  in  tha  various  casas  of  cclloca- 
■+■400  oxc  « ^4 in  +’«.^ it.  • 
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Ca  s«  i 


?a  ta 


I  Jnraducao. 

ir  ^aducad  aith  topographic 

:^ata  sat  C  and  Rapp  160. 


Ill 


^aducad  jith  Tocograohic 
data  sat  ^  and  Rapo  180. 


:v,v 


^a.*«  as  III 


VI 


Raducad  only  aith  topographic 
data  sat  F. 


Table  6.  Identification  of  collocation  cases. 


Tha  rasults  o*  collocatLon  ara  shoan  in  tabla  7.  Wa 
can  saa  that  tha  rasplts  for  jrawity  anomalias  ara  signifi¬ 
cantly  i.-aprovad  by  usin;  tha  obsarvations  raducad  aith  son* 
topographic  data.  7na  usa  of  tha  spharical  harmonic  coaffl- 
cients  doas  not  affact  tha  rasults  for  gravity  anomalias. 
Tna  bast  laval  of  accuracy  obtainjd  in  this  t«st  aas  3  mil- 
ligals  (ora  sigma).  This  rasuLt  »ss  obtainad  oy  using  300 
gravity  anomaly  cosarvati  ons«  soecad  about  5  .ninutas  apart, 
fha  inclusion  of  oaflaction  or  unculation  obsarvations  doas 
not  appaar  to  affact  tha  rasults  ♦or  gravity  anomalias. 
*^aanufhila«  by  incraesing  tha  numoar  of  ancmaly  obsarvations 
>y  a  factor  of  3»  aa  gam  apout  ona  milligal  in  accuracy. 

For  daflactionst  tha  bast  ?ccuracias»  again  in  tarms  of 
standard  eaviation  of  tha  arrors*  for£  u.»vi  r|  o-rt  o.hq  arc. 
^ec-O'hds  0.^9  arc. 
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obsarvations 


(sbiervation  -  pradlctad) 


Casa 

n 

1)8  X 

isa  an 

std.  dav 

1 

A9 

100 

400 

20.25 

-1.21 

6.30 

3 

20 

2.59 

-3.62 

1.11 

*? 

3 

20 

-3.13 

-0.61 

1.22 

N 

0 

6 

-27.35 

-25.00 

0.50 

II 

AS 

100 

400 

-15.  37 

-0.14 

4.13 

i 

3 

20 

-3.96 

-0.12 

1.34 

V 

3 

20 

5.21 

3.25 

1.54 

N 

0 

6 

-  J.25 

-3.00 

0.50 

III 

iOO 

400 

-16.37 

-0.14 

4.15 

3 

20 

-3.44 

0.12 

1.36 

n 

3 

20 

2.68 

0.65 

1.13 

IV 

As 

300 

200 

-12.52 

0.26 

3.00 

f. 

20 

-3.81 

0.16 

1.39 

n 

3 

20 

1.97 

0.12 

0.97 

V 

As 

300 

0 

5 

1 8 

2.s7 

1.13 

0.92 

V 

18 

-2.01 

-0.08 

0.94 

v: 

As 

3C0 

2  00 

-12.51 

0.23 

3. 03 

? 

3 

20 

-4.91 

-1.51 

0.79 

1 

3 

20 

-l.<»9 

0.09 

0.8  9 

Table  7.  Collocation  results.. 


graphic  and  isostatic  radLCtions,  naglacting  th»  Rap?  160 
conponants.  Tht  maximum  and  naan  arrors  for  th#  north-south 
conponant  ara  sli_shtly  grsatar  hoaavsrt  than  an«n  th#  <5pp 
180  coaffieiants  ara  usad.  Tn#  banafit  of  using  tha  fiva 
sacond  innar  grid  data  is  liKanisa  to  oring  tha  aaan  arrors 
closar  to  zaro*  tha  standard  daviation  ranaining  tha  sama. 

Tha  casa  of  undulations  is  avan  naraar  to  arslyza« 
givan  only  six  ob sarvations.  From  this  tast*  aa  aould  con- 
cluda  that  tha  bast  aay  to  astimata  thasa  is  strictly  from 
tha  spharical  harmonic  :o af f xci »n ts>  Tha  a^facts  of  raduc- 
ing  tha  obsarvatlons  and  applying  collocation  dagradad  tha 
rasults  obtainad  only  aith  tha  soharical  harmonics* 
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CONCLUSIONS  AND  RECOMMENDATIONS 


A  t»chniou»  ^or  computins  a  local  s»*ewity  fiald  mocal 
in  mountainous  tsrrain  has  baan  tastad.  Tha  rasults  indi- 
cata  that  gravity  anomalias  can  oa  cradictad  aitn  an  accu¬ 
racy  of  3  millijals*  daflactions  of  tha  vertical  with  an 
accuracy  of  about  1  arc  second  and  gaoid  undulations  aith  an 
accuracy  of  about  1.0  maters. 

For  gravity  anomaliast  la?st  squares  collocation  fol- 
louing  trend  removal  by  topo gr apnic/is os t a t i c  data  and 
spherical  harmonic  coefficients  gave  the  bast  results.  For 
deflections*  tha  inclusion  of  tha  spherical  harmonic  ccm- 
nonant  haloed  reduce  ?  bias  in  the  errors.  For  undulations* 
the  bast  results  were  oCtairec  directly  from  the  spherical 
harmonic  coefficients*  neglecting  any  tocoshephic  data  or 
least  sauares  collocation. 

Tna  bast  comtainaticn  of  terrain  data  acoesrs  to  be  a 
lidely  spaced  (eoout  I  degree)  outer  grid  extendinr  to  ?t 
least  10  degrees  around  the  test  points  ena  ct  least  one 
coarser  inner  grid  Caoout  5  minutes)  extending  to  eoout  5 
dagraas  from  tha  computation  point.  The  rasults  are  signi¬ 
ficantly  annancac  oy  interpolating  tnie  inner  grid  to  an 
even  finer  grid  in  the  vary  near  vicinity  of  the  computation 
ocmt.  Wnen  this  tyoe  of  interpolation  is  cone*  it  may  ba 
unnecessary  to  include  a  very  fine  grid  of  elavations  as  a 
third  "inner-inn er"  zcna. 
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Th*r«  sKist  v«ry  afficient  computar  projranis  for  ?ei — 
forming  sll  of  th«  calculations  dascribsd  narcin.  Givan  a 
madium  siiad  irinicomputar  i  ail  of  th#  calculations  can  ba 
carriad  out  in  a  faa  hours.  Tha  tocogr ^pnic/ isos tatic  com¬ 
putation  is  tha  most  t im a-cons um in g  out*  using  appreximate 
formulast  thasa  can  ba  dona  in  a  raasonabla  time. 

For  adjusting  triangulation  nataorks*  or  reducing  hor¬ 
izontal  and  vertical  angles  to  tha  ellipsoid*  the  accuracy 
of  tha  vertical  deflections  obtained  herein  is  sufficient 
for  even  first  ordar  work.  If  tha  elevation  is  7  degrees* 
an  error  in  deflection  of  tha  vertical  of  2  arc  seconds  aill 
result  in  an  error  of  0.25  seconds  in  norizcntal  angles* 
ahich  is  beloa  the  observation  arror  o*  most  angles. 

In  mountainous  areas*  conversion  of  satellite-derived 
neights  to  heights  above  sea  level  reguiras  tne  geoid  undu¬ 
lation.  This  study  indicates  that  these  can  ba  computed  to 
an  accuracy  of  about  1  metar*  ahich  is  consistent  aitn  pest 
results.  This  will  propogate  directly  into  a  1  meter  error 
in  sea-level  heights*  naturally.  In  tha  absence  of  other 
information*  this  may  be  an  acceptable  error  far  some  jiork . 

The  accuracy  of  the  gravity  anomalies*  3  milligals* 
means  that  they  coula  be  used  »for  instance,  to  compute 
deflections  of  the  varticsi  and  gaoid  undulations  through 
Vening  '^einesz'  and  Stakes'  integrals  respectively.  The 
predicted  anomalies  coula  also  be  useful  in  geoohysical 


axplorations 


For  alignmant  and  corrictioo  af  martial  navigation 
systana*  tne  accuracias  ootainad  for  daflactions  and 
anomalias  could  oa  uaaful  in  certain  applications.  This 
araa  must  ba  exclorad  ♦urthar  nomavar. 
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